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Abstract 

Classical Edgeworth expansions provide asymptotic correction terms to the Central Limit 
Theorem (CLT) up to an order that depends on the number of moments available. In this 
paper, we provide subsequent correction terms beyond those given by a standard Edgeworth 
expansion in the general case of regularly varying distributions with diverging moments (beyond 
the second). The subsequent terms can be expressed in a simple closed form in terms of certain 
special functions (Dawson's integral and parabolic cylinder functions), and there are qualitative 
differences depending on whether the number of moments available is even, odd or not an 
integer, and whether the distributions are symmetric or not. If the increments have an even 
number of moments, then additional logarithmic corrections must also be incorporated in the 
expansion parameter. An interesting feature of our correction terms for the CLT is that they 
become dominant outside the central region and blend naturally with known large-deviation 
asymptotics when these are applied formally to the spatial scales of the CLT. 

1 Introduction 

Let (Xfc : k > 1) be a sequence of independent and identically distributed (iid) random variables 
(rv's) with mean zero and finite variance cr^. If Sn = Xi + ... + Xn, then the standard Central Limit 
Theorem (CLT) says that for each fixed a; G M, 



P (Sn/n^^^ <x) =P{N (0, a^) < x) + a (1) 



as n oo. Edgeworth expansions provide additional information about the previous error term, 
but only if higher moments exist. In particular, if E\Xk\"^ < oo for m > 3 and assuming that 
the Xfi-^s are continuous random variables (or simply non-lattice) an Edgeworth expansion of order 
[m — 2j takes the form 

[m-2j 

P [Sn/n'/^ <x)=P{N {0,a^) < x) + ^ n'^/^p. ^ ^ 2J/2 

as n 00, where Pj (•) is appropriately defined in terms of Hermite polynomials and involves 
coefficients that depend on EXj^"^ , for fixed a; G M (see Section 2 below). As a consequence, 
a traditional Edgeworth expansion cannot continue beyond the [m — 2j-th correction term in a 
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situation where E j = oo. Our focus here is precisely on this situation, building on 

the preliminary results of Bazant (2006). 

In particular, our goal in this paper is to provide additional correction terms to the Edgeworth 
expansion and further reduce the order of the error term as n oo in the presence of regularly 
varying distributions (which are basically power laws; see Equation ([2])). As we shall discuss in detail 
below, an interesting feature of our development is the fact that the additional correction terms 
are qualitatively different (giving rise to expressions involving log (n) in some cases) depending on 
whether the power-law decaying rate of the distribution is an integer number or not, or, when such 
rate is an integer, on its parity. 

After this paper was completed we learned that a related set of expansion is given in Chapter 2 of 
Vinogradov (1994). Our results here appear to be more general than those expansions. In particular, 
our results apply to regularly varying tail of any real number a > 2, while Vinogradov considers 
only power-law tail with integer a. Also, our analysis allows us to obtain explicit expressions for 
the coefficients in our Edgeworth expansion in terms of special functions, an important interest in 
this paper that was not pursued by Vinogradov. 

Our results can be applied in the context of random walks (Hughes (1995); Weiss (1994)), 
where corrections to the CLT are sometimes called Gram-Charlier expansions (and extended to 
vector random variables). As is typical for such expansions, our asymptotic formulae can be very 
accurate for small n, but unlike standard Gram-Charlier expansions, our corrections also become 
dominant at the edge of the central region and beyond, as illustrated by numerical examples in 
Bazant (2006) (see also Theorem [3] in Section 4). Random walks with heavy-tailed distributions (or 
synonymously "fat-tailed distributions" in physics and other literatures), in particular those that 
are regularly varying, are increasingly used to describe stochastic processes in diverse fields, ranging 
from physics (Bouchaud and Georges (1990); Metzler and Klafter (2000)) to economics (Bouchaud 
and Potters (2000); Embrechts, Kluppelberg and Mikosch (1997)). 

Edgeworth expansions are also often used in statistics and applied probability. For instance, 
in statistics these types of expansions have enjoyed successful applications in statistical analysis 
involving relatively small samples (Field and Ronchetti (1990)). In bootstrapping, where data are 
resampled to calculate quantities related to the estimator, Edgeworth expansion can reduce the 
order of the approximation error (Hall (1992)). In applied probability settings such as queueing or 
insurance one is often interested in approximating performance measures of a stochastic system, for 
instance the tail of the delay in a queueing system or the ruin probability in insurance, by means 
of a CLT. In some cases, one can derive corrections to such approximations (Blanchet and Glynn 
(2006)). These corrections are often called corrected diffusion approximations and are similar in 
spirit to Edgeworth expansions. In fact, our ideas here can be applied to developing higher order 
corrections to a class of queueing models in heavy-traffic; such development will be explored in 
future work. 

Other types of approximation that are often applied to the distribution of 5„ include so-called 
large deviations or moderate deviations asymptotics. In contrast to Edgeworth expansions, which 
approximate the distribution of Sn in spatial scales of order O (n^/^) , large deviations results typ- 
ically provide approximations in spatial scales of order O (n), whereas moderate deviations are 
derived in spatial scales of order O (n^/'^^^) for e S (0,1/2). In situations where the increment 
distribution has finite logarithmic moment generating function, standard results provide an expo- 
nential decay rate to the tail of Sn in large deviations or moderate deviations scales, say for quantiles 
that are of the form nx. Using the fact that typically the rate function is twice continuously differ- 
entiable at the origin (as a function of x), one can relate (or transition from) large, moderate and 
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CLT approximations smoothly. On the other hand, when the increment distributions are regularly 
varying, the large deviations rate, as we shall recall in future sections, turns out to be basically 
polynomial. Given that further terms in an Edgeworth expansion reduce the polynomial (in powers 
of v}/'^) error rate of the tail estimate, one may wonder how Edgeworth expansions relate to large 
deviations rates as additional correction terms are added. As we shall see, perhaps surprisingly, 
our additional terms blend smoothly with large and moderate deviations rates when these rates are 
formally applied to spatial scales of order n^/^ only when the power of the regularly varying tail 
is non-even. These formal connections to large deviations asymptotics are further explored in Sec- 
tion 4; there we would discuss our asymptotics in relation to the classical uniform large deviations 
result by Rozovskii (1989), in which the sum of iid regularly varying rv can be approximated by a 
sum of Gaussian distribution and heavy-tail asymptotic. We will also show that under restrictive 
conditions our asymptotic provides a refinement of this result in the moderate deviations region. 

As we indicated earlier, our development gives rise to qualitative differences in the form of the 
correction terms depending on the decay rate at infinity of the increment distributions. In order 
to further explain such qualitative differences, let us introduce some notation. We impose the 
assumption that the XkS have a positive density fx (•) of the form 



X [X) 



(1 + |x|)-("+^)l(|x|), 



where L (•) is a slowly varying function (i.e. L (xt) /L [x) — > 1 as x oo for every t > 0) so that 
var{Xj^) = 1 (our development also deals with the non-symmetric case; see Section 2.2 and 5). Also 
assume that \(j){9)\d9 < oo where (j){9) = Ee^^^*" is the characteristic function. Denote t/(x) as 
the standard normal density. Theorem [2] in Section 2 yields that the density of Sn/n^^"^ behaves 
asymptotically as 



1 + 



Gj{x) 



3<j<a 
j even 



L(n^/^x) 



+ 



'^G^{x)+o{ 

Gaix) + O 



„a/2-l ) 
\ „a/2-l / 



when a is not even 
when a is even 



as n oo, where Cl(^) ^ suitable slowly varying function (see Equation ([6])) and Ga{x) is written 
in terms of special functions whose qualitative behaviors depend on the parity of a (see Equation 

m)- 

Our approach to developing the Edgeworth expansions is standard. First we develop an ap- 
propriate Gram-Charlier series for the Fourier transform and then apply inversion. An interesting 
element of the analysis in our situation, however, involves the contribution of a non-analytic com- 
ponent in the Fourier transform. We shall build the elements of our expansion starting from the 
symmetric case and then extend the construction of our approximation to more general regularly 
varying densities. Finally, we note that our approach can be used to develop asymptotic expansion 
for stable limits, which will be reported in future research. 

The paper is organized as follows. First, we discuss our main results in Section 2. Then we 
discuss in detail the symmetric case in Section 3. The formal connection to large deviations is given 
in Section 4. Section 5 studies the case of non-symmetric density. 
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2 Main Results 

2.1 Review of Ordinary Edgeworth Expansion 

Let (Xk : A; > 1) be a sequence of iid rv's with mean zero, unit variance and with a regularly varying 
density. More precisely, we assume that the density / (•) of Xf^ satisfies 

/(x) = {l+xy^^^^H+{x)I{x>0) 

+ (1-x)"(^+^)l_(-x)/(z<0), (2) 

where /3,7 > 2 and (•), L_ (•) are slowly varying functions at infinity (i.e. L± (cb) /L± (b) — > 1 
as b — > oo) defined on M. We shall write <j) (9) = Eexp {iOXj.) for the characteristic function of X^ 
and put il) {0) = log (j) {&)■ Assume \(j){9)\d6 < oo. Finally, we let Sn = Xi + ... + X„, n > 1, be 
the random walk generated by the X^s, set {9) = E exp [OSn/n^^'^) = 0" and define 

V', (0)=log0„ (0)=nV^ (0/nV2). 

As we indicated in the introduction, our goal is to extend the standard Edgeworth expansion 
given by approximation ([T]), which is obtained from a formal expansion of the Fourier transform of 
Sn/n^^"^- In particular, since V'' (0) = we have, for fixed ^ G M, 



, (n, i/2\ ^ ^ji^oy / \ 



j=2 

if m < min (/3,7), where Kj's are the cumulants of X^. Therefore, 

= exp(-«V2)»p(f:Mg. + „ 



( 



exp (-^V2) 



^ ^ £ fe! ( £ 7!nJ/2-i 




( ra q / 1 " 

\ 5=3 fc=l ^ 



(3) 



where is the coefficient of in the expansion of (Y^'j^^ f^jiidy /{jW^"^ "'^))'^- Hence 

the summation in ([3]) has first term K3 (i^)'^ /n^/^^ second term (^3(20)^/72 + K4(i0)^/4!)(l/n) and 



so on. 



Recall that the inverse Fourier transform of (—iO) exp (—0^/2) equals d^r] (x) /dx^ , where 
r] (x) = exp (— x^/2) /\/27r. In turn, we have that 

d^ 

-^rj {x) = 7] (x) (-1)'' Hk (x) , 

where (x) is the A:-th Hermite polynomial. In particular, Hi (x) = x, H2 (x) = x^ — 1 and 
H3 (x) = x^ — 3x etc. Therefore, taking inverse Fourier transforms in both sides of ^ we obtain 
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(see Hall (1992) for a derivation of the distribution function, and Evans and Swartz (2000) and 
Section 3 for density function) that the density (•) of Sn/'n}^'^ satisfies 

m 

fsjn^/- (^) = v{x) + ri {x) J2 n-'/'+'Gk (x) + o (n"™ Z^+i) , (4) 

fc=3 

where Gk (•) is a polynomial of degree k that depends only on the first k cumulants of Xj. In 
particular, we have that G3 = K^Hs/Sl, G4, = kIHq/72 + ^4^^4/4!. 

Our goal is to understand the contribution of the error term o (n~™/^+^) in expansion (jH and 
the problem is that subsequent G^'s (for k > min(/3,7)) involve moments (or cumulants) that are 
not defined for Xj. The idea is to write ip (0) = x (d) +C (^) where x (■) is analytic in a neighborhood 
of the origin while ^ (•) is a non-analytic function. Dealing with the analytic component, x, gives 
rise to the standard Edgeworth expansion as in so the interesting part involves the analysis of 
^ (•) which, as we shall see, yields Theorem [J and [2] that extends 



2.2 Main Theorems 

Throughout this paper we assume the following differentiation condition on L± (for simplicity we 
write L below): 



Assumption 1. L is bounded, dijferentiable and satisfies 

'L{x) 



L'{x) = o(^- 



(5) 



as X ^ 00 



Assumption [T] controls the order of the derivative of L. Examples of slowly varying function 
that satisfy ^ include logx, (logx)^ etc. 

We also introduce the function C,l{x) defined as 

c,{x)=r^du (6) 

Ji u 

Note that Assumption [1] implies 

L{x) = o{Cd^)) (7) 

as X — )• 00. Both Assumption [1] and this condition will be used in our proof of the main result. 
Next let us define Dawson's integral D[z) as 

D{z) = e-"' r e'^dt (8) 
Jo 

and the classical parabolic cylinder function Dy[z) with parameter where Re(z^) > —1, as 

D,{z) = ./^e^'/4 r e-*'/2t- cos{zt - —)dt (9) 
V vr 7o 2 

We state our main theorem: 
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Theorem 1. Assume (X^ : k > 1) are iid rv's with zero mean, unit variance and density following 
^ with L-i- and L_ satisfying AssumptionUl Also assume \<j){9)\d9 < oo. We have, for fixed 
X > 0, 

fs^/n^ni.^) = ri{x) 

as n ^ CO, where Gj{x),j < min(/3,7) are the ordinary Edgeworth coefficients (as discussed in 
Section 2.1) and F{x,n) is defined according to the following cases: 
Case 1: /3 = 7 = q 



1 + 



Gjix) 



^ ^i/2~l 
3<j<min(/3,7) 



+ F{x, n) + o{F{x, n)) 



F{x, n) 



r(Q+l)n°'/2-l 



H^{x) [CL^{ny^x)+CLM'^^)) 



'V2 



r(a+l)n°'/2-l 



^^£^D{^){L^{ny^x)+L_{ny^x)) 



■V2 



2 r(o+l) sm(a7r) 

+ Da{-x)L_{n^/'^x 



-^^[D^{x)L+{ny^x) 



for even a 



for odd a 



for non-integer a 



Case 2: (3 <-f 



Fix, n) 



Case 3: (3 > J 



2 r(/3+l) sm(/37r) n/3/2-i 



Fix, n) 



2^r(7+i) nT/2-1 

2^r(7+l) nT/2-1 
2 r(7+l) sm(77r) 



for integer fi 
for non-integer (3 

for even 7 

for odd 7 

for non-integer 7 



Here Hk{z), D{z) and D^^z) are Hermite polynomial of order k, Dawson's integral and classical 
parabolic cylinder function with parameter v respectively. 



Note that this theorem states the result for x > 0. For x < 0, one merely has to consider — 
and — 5n/n^/^ and the result reduces to the x > case. 

This result can be obtained by first considering a symmetric density and then splitting the 
non-symmetric density into odd and even functions which can be tackled separately using result for 
the symmetric case. For the symmetric case, i.e. /3 = 7 = a and L = = L_, we have a neater 
representation: 
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Theorem 2. Assume (X^ : A: > 1) are iid rv's with unit variance and symmetric density f{x) = (1+ 
|x|)~(""'"^''L(|x|) with L satisfying AssumptionUl and a > 2. Moreover, assume \(l){9)\d9 < oo. 
Then, for fixed x £M, 



/s„/nl/2(2;) =r/(x) 



„i/2-i 

3<j<a 
j even 



+ 



L{r^!2llr (r) I nf ^("'^'^) ^ 



,a/2-l 



for non-even a 
for even a 



as n oo, where Qj^ is defined in Here Gj{x),j < a are the ordinary Edgeworth coefficients 
(as discussed in Section 2.1), while Ga{x) is defined as 



V2 d° 
' r(a+l) dx' 



T2> 



for odd a 



[Da{x) + Da{—x)] for non-integer a 
iri^^Haix) for even a 



2 r(p+l) sin(a7r) 

2 e— /2 



(10) 



where D{z), Dy{z) and H^^z) are Dawson's integral, classical parabolic cylinder function with pa- 
rameter V and Hermite polynomial of order k respectively. 



As we indicated before, the analysis behind the previous result involves understanding the 
behavior of the non-analytic component ^. In fact, we obtain the following decomposition of ip: 

Proposition 1. Assume a symmetric density as in Theorem\^ We have 

m = x{.o)+m+o{m) (n) 

where 

l<j<a •'' 
j even 

is the ordinary Taylor series expansion up to the largest moment and 

{— TV — rn-^ — — I^I^^Crar) for non-even a 
r(a+i) W^L{\e\) for even a 

is the non-analytic component. Here is defined in 



With this decomposition, we can perform similar derivation as Equation ([3]) to obtain 
i?exp(^5„/nV2) = exp (-0^2) ( 1 + f E ^ W^);^ + < (;^) + « (< (;^) ) ) ^ 

\ Q — 3 k — 1 / 

Dawson's integral and the other special functions then arise as the inverse Fourier transform of 
g-s /2|^|a appearing in the non-analytic term of the above expression. 

The rest of the paper provides the necessary details behind the analysis of ^ and the Fourier 
inversion required to obtain Theorem [2] and then Theorem [U the proof of Proposition [1] follows 
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directly from Proposition [3l which is given in Section 3. As we indicated in the introduction, in 
Section 4 we shall discuss the formal connection between our expansion and that of Rozovskii (1989) 
in the context of large deviations. It turns out that the extra term in our expansion, which comes 
from the non-analytic component of the characteristic function, blends smoothly with the heavy-tail 
asymptotic, at least in some restrictive cases. When this holds, there exists a smooth transition 
from central limit to large deviations region. 



3 Symmetric Density 

Our goal in this section is to prove Theorem [2j As noted before, we adopt the approach of first 
expanding the characteristic function of a symmetric regularly varying density, and hence its cu- 
mulant generating function, followed by an inverse Fourier transform to get the stated result. In 
order to obtain the characteristic function for regularly varying density, it helps to first consider in 
detail the nature of the approximation in a simplified context, namely, that of a symmetric Pareto 
density. 



3.1 Characteristic Function 
3.1.1 Pareto Density 

In this section we shall assume that 

where a/ is the normalizing constant that makes / (x) a well defined density. For example, a/ = 
7r/\/2 if a = 3. With this we have the simplest version of our expansion: 

Proposition 2. For rv following {14% the characteristic function takes the form 

^(«, = 1 + V t}y^e' + 1 -ff^ ""° ^ " (15) 

25?„ J- I - r('!+l)' |9riog|»|+o(|g|°log|D|) forevena 

j even 

as —7- 0, where mj is the j-th moment of the rv. 

The proof of this theorem is divided into four cases: a is odd, even, and non-integer with the 
integral part being odd and even. Each case involves finding the asymptotic of an integral as below. 
Also note that since the density is symmetric, the characteristic function is real and symmetric. 
Thus we can consider only 9 > without loss of generality. Moreover, throughout the paper we shall 
use the notation "~" to denote equivalence relation i.e. f{x) ~ g{x) iff f{x)/g{x) — )■ 1 as x \ (or 
X oo depending on context). We also use \x\ as the integral part of x i.e. max{A; G N : A: < x} 
and [x] as the ceiling of x i.e. min{A; S N : /c > x}. 

Lemma 1. For ^ \ 0, we have the four asymptotics: 
Case 1: a is odd 

^ 'dxr^-ire (16) 



1 + X 



l+a 



Case 2: a is even 



foo ^a—l„iOx 



_|_ |7-|1+" 

For non-integer a, denote q = lqj. Then 
Case 3: q is even 



/oo ^ . 
r^—dx^-2ieioge 
1 I I ™n -4~rv ^ 
-OO 



Case 4' Q is odd 



Proof. We prove the lemma case by case: 

• Case 1: a is odd 

Since a is odd, the integrand is a symmetric function, and we can write 

/oo a—l/ idx _ -[\ fOO ™a— 1 
^, ' dx = 2 / — ^— _(cos(0x) - \)dx 

roo u°'~^ 

= 20 I — rr ; icosu — l)du 

by a change of variable it = 0a; in the second equality. Since 

(cos u—1) 



< ^1 cos n — 11 



which is integrable, by dominated convergence theorem we get 



■"-^ '•^COSW-l, TT 



I 7"; — (cosu — — )• / 



-du 



^2 2 

as ^ 0. Hence Case 1 is proved. 

• Case 2: a is even 

Note that the integrand is an odd function, and 

r-oo ^a-i^Wx /•°° x"-^ sin(6'x) 



/oo ^a—L^itlx poo ^a- 
r-rn-dx = 2i — 
1 + Jo 1 



+ x 



1+a 



-dx 



2i9 / ; — du 
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again by a change of variable u = 9x. Note that as — )• 0, the integral blows up at n = 0. So 
for small e, we write 

/"°° u'^-^sinn , f n^-^sinu , u'^'^smu , 
; — du = / ; du + / au 



Now 



du < I —du < oo 



On the other hand, 

1 



/ -TT ^ — du = / -TT ^ — + i?(6') 



where 

,2+a 



1^(^)1 < C / ; — du <C udu = 0(l) 

' ^ ^' - Jo + - Jo ^ ' 

for positive constants C and 

r 1^ , t^u = log(^i+" + «^+") ' - - log e 

Hence Case 2 is proved. 

• Case 3 and 4: a is non-intcgcr with q even and odd 

The proof is similar to the last two cases. For even q, we write 

= 20°^-« / ^- ^ — (cos-u-lWu 

'/ 

JO 



and for odd q, we write 



/oo x'^e x'' 
1 + Jo 1 



1+a 



+ X 

oo „,q„- 



dx 



sm u 



= 2i^"~^ / ^- : du 

Jo ^1+" + 

Jo 

by using a change of variable u = 9x and dominated convergence theorem. We use the fact 
that 

/•oo 

/ u^-^(e^" - l)du = T{z)e'^ for - 1< Re(z) < 

to get 

f°° cosu — 1 , ^, , ,— a + (7 , 

/ — 7- au = r(— a + o) cos( — tt) 

io u^+<^-i ^ ^' 2 ' 
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and 



/■°° sinu -a + q 
/ — n an = 1 —Q + g sm vr) 



The lemma is then immediate. 

□ 

We are now ready to prove Proposition [21 The proof basically reduces psp into expressions in 
terms of the above four integrals through multiple use of L' Hospital's rule. 

Proof of Proposition\^ It suffices to consider > 0. Again we consider case by case: 

• Case 1: a is odd 

Applying L'Hospital's rule successively for a — 1 times, we have 

POO ate'-^'^ J 1 (-'^y^^i^i ai 

J-oo T+t^^^ - 1 - El<i<« j! 
T j even 

1™ T^S 

J even 

= lim 

Loo - E2<i<a (^2)! ^ 

J even 



lim 



a-2 



>o a(a - 1)6' 



Irni ^-^ 

_1 ("-i)/2o ' dx 
lim 



a! 0-s>o 

(_l)(»+l)/2^^^ 

a! 

where the last step follows from (|16p . Note that the integral in each step above is finite since 
the moments exist up to order a — 1. 

The result for odd a thus follows. 

Case 2: a is even 
Using the fact that 

^fe" log A; = a(a - 1) ••• (q - m + 1)A;— log k + Ta-^A:""™ 
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for m = 1, . . . , a — 1, where Ja-m constants, and applying L' Hospital's rule for a — 1 
times, we get 

POO a^e'^^ , -I {-'^V ^'^'^i nj 

J-oo l+\x\^+-' dx - 1 - El<i<a j! 

1 even 

lim ^ 

e^o riog^ 

TOO g/i^-^z^-^e''^'' ^ 

— iim 



9^0 a\9\og6 + ^i6 



a/z° J-oo T+CTn^"^ 

= — ; — lim — 

a\ 6 log 6 

2(-l)°/2a/ 

a\ 

where the last equality follows from ([IT]). 
The result for even a thus follows. 

Case 3 and 4: a is non-integer with q even and odd 

Following the same line of proof as the previous two cases, we get, for even 

J-oo TTN^^^ - 1 - Ei<i<" i! '' ^^ 

1- 7 even 

1™ TJS 

■g roo xg(e'9"-l) , 
* "/J-oo "-^ 



lim 

9^0 a(a - 1) • • • (a - g + 1)^"'' 

(-l)'?/^a/r(a-g + l) + g . 

; 2r -a + g COS — vr) 

T{a + 1) 2 

2(-l)''/27ra/ cos(^^) 



r(a + l)sin((-a + Q')7r) 

T{a + 1) sin(^^7r) 

r(a + 1) sin ^ cos f 
r(a + l)sinif 



and for odd q, 



13 



TOO aje^^'' ^ 1 (~l)^^^™i/3i 

J-oo l+U.|l+""^ - 1 - 2^1<j<Q j! 

1. 7 even 

1™ TiS 



lim 



g POO ^r^e!!^j^ 
^^a/ J_oo 



>Oa(a- l)---(a-g + l)r-« 

(_l)to+l)/2 r(a-g+l) , -« + g . 

• 2T[—a + q) sm( vr) 

r(Q + l) ^ ^' ^ 2 ' 

2(-l)(9+i)/2^a/sm(^7r) 

r(a + 1) sin((— a + q)T:) 

(_l)(9+l)/2^a^ 



r(Q + 1) cos(^^^7r) 

(_l)(g+i)/2^a/ 
r(a + 1) sin sin f 

r(a + l)sin^f 



by using the reflection property of gamma function and also the double-angle and sum-of-angle 
trigonometric identities. 

We have proven the proposition. 

□ 

3.1.2 Regularly Varying Density 

We are going to extend the result in Proposition [2] to density with regularly varying tail in the form 

/(^) = rT]^' "^^^""^'"^^ ^^^^ 

where L is a slowly regularly varying function and satisfies the differentiation condition in Assump- 
tion [H 

Proposition 3. For rv following hl^) with L satisfying Assumption\^ the characteristic function 
takes the form 

6(9) = 1 + V tD^QJ + / ~ r(.+ilsin^ l^r^ (w) + (l^l"^ (w)) ^"'^ ^"^-^^^^ " 
2ir<a \ ^(S? (^) + o {m, (^)) for even a 

j even 

(19) 

as ^ — 7- 0, where mj is the j-th moment of the rv and Cl is defined in 



To prove Proposition [3l we need the following property of slowly varying function: 
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Lemma 2. Let L be a slowly varying function. Then for any Pi, p2 > 0, we can pick an r] > such 
that for0<6<r], 

L{§) ^ r CxP^ for x>l 
7^(1) - \ Cx-P2 for X <1 

where C is a positive constant. 



Proof. Karamata's representation (see for example, Resnick, 1987) states that for a slowly varying 
function L : — t- M4., we have 



L{x) = c{x) exp 1^ t~^e{t)dt^ 



where c : M+ — t- M+ and e : M — t- M such that 

lim c{x) = c E (0,00) 



lim e{t) = 



So we have 



Since L is bounded, we can choose c{x) bounded from above and c{l/9) > ci > for 6 small 
enough. So 

c(J) - ' 

for some C > for small 9. Also by a change of variable t = s/6, we have 

exp jy^" t"^e(t)(it| = exp |^ s~^4^)ds^ 

Suppose X > 1. For any given pi, if we pick 6 small enough, we have e(|) < pi for all s > 1, and 

exp I s~^e(-)(is| < exp = 

The result for x > 1 follows. Similarly, for x < 1, if we pick 9 small enough we have e(|) > —p2, 
and 



exp 



s ^e(^)(is| = expj-y^ s ^e{^)ds^ < exp I^J s ^P2ds 



Thus the lemma is proved. □ 



We now state the correspondence of Lemma [T] for regularly varying density: 
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Lemma 3. For 9 \ 0, we have the four asymptotics: 
Case 1: a is odd 



°° x"-iL(|x|)(e^^^ - 1) , fl 
-ax ~ —ttBL — 



l + \x 



l + OL 



Case 2: a is even 

foo ^a—lTn„\\AOx 



where is defined in For non-integer a, denote q = loj as the integral part of a. Then 



Case 3: q is even 

„i6x 



Case 4- Q is odd 



Proof. Again we prove the lemma case by case. 
• Case 1: Q is odd 

We follow similar line of proof as in lemma [TJ Using a change of variables u = Ox we obtain 

= 261 / -r— — 7— -(cos n - Ijdu 



■oo 



l + |x|i+° io 0^+° + ni+^ 

.1. /■'^ n°-^ 



2^^(77) / ■ (cosu-l)d 







Now using lemma O for small enough, we can pick < pi,/92 < 1 such that 



-J- • (cos n — 1) 



C C 
< „ , (cos n - l)/(n < 1) + ^5 (cosn - l)I(u > 1) 

^/+P2 Pi 

for some C > 0, for small enough 9. Note that the majorizing function is integrable. Also 
note that by the definition of slow variation we have 

as — )■ for all u > 0. So by dominated convergence theorem, 

\ — n ; — • — TT- • (cos u — llau — )• / ^(cosn — llan 

and the asymptotic follows. 
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• Case 2: a is even 
We write 



1 + |x|i+" 



dx = 2i 
= 2ie 
= 2ie 



oo 



L{x) s\n{9x) 







1 + 



./o 



.a— 1 r (' " ^ 



dn + 



oo yOi—l T (u 



sinti 



9 



l+OL 



+ u 



1+a 



du 



for some small e > 0. For the second integral, if we choose < < 1 in Lemma [21 we have 
the following bound 



~ u"-iL(^)sinu 



-du 



for some C, Mi > 0, for small enough 0. 
For the first integral, we write 

LL(^)sinn 



<L|i 



1+a 



du < MiL I ^ 



(20) 



1+a 



du 



where 



\R{e)\<M2L(^ 



^l+a _^ ^1+a 
^a+2^-p2 







du + R{e) 



^l+a ^ ^l+a 



du < KhL - 



(21) 



for some M2, M3 > 0, < P2 < 2 and small enough. 
Now do a change of variable x = u/6 back, we get 



gl+a 



,1+a 



-du 



x"L(x) 
1 + xi+" 



dx 



L(x) 



dx 



(22) 



Combining (j20p . (j2ip and ()22p . we obtain the result for even a. 

• Case 3 and 4: a is non-integer with q even and odd 

Using the same line of proof as Lemma [1] and Case 1, and taking < < a — q and 
< P2 < 2 — a + we get the stated result. 

□ 

The proof of Proposition [3] is similar to that of Pareto density, and hence we only outline the 
proof below. Proposition [1] then follows immediately by letting ij: = log (p and applying a Taylor 
expansion. In doing so one just has to verify that the non-analytic term in (f) is the only term of 
that order in the expansion of ip to conclude the result. 
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Outline of Proof of Proposition\^ Again it suffices to consider > 0, and we use L' Hospital's rule 
successively. For odd a, recognizing that L'{^) = o{0L{^)) as — )• by Assumption [H we have 

POO L(x)e^»- , , ^ (-ly/^m. ^j 

J-oo ~ ~ Z^l<J<a j! 

V ) even 

iim 



lim 



lim 





rL(i) 






TOO 

J—oo 


j even 














roo 
J—oo 


J even 




QJ-l 



_1 («-l)/2 -'^^ 

= ; iim -i- 

a! 9^0 6IL(i) 

(_l)(a+l)/2^ 

The case for even a can be proved similarly by using (HD that states that = o(Cl(^)), while 
the case for non-integer a is proved similarly with Assumption [TJ 

□ 

3.2 Beyond Edgeworth Expansion 

We will now prove the main theorem for the symmetric case, namely Theorem [2j First of all we 
need to evaluate the following integrals, which corresponds to a Fourier inversion of the non-analytic 
component of the expansion studied in Section 3.1: 



Lemma 4. 



oo 



' 2^/2(-l)(— for odd a 
Y^sec ^e~^^/'^[Da{x) + Da{—x)] for non-in. 
{—\)^l'^\p2Me~^'^l'^lia[x) for even a 



where D{z), D^{z) and Hj.[z) are Dawson's integral, classical parabolic cylinder function with pa- 
rameter V and Hermite polynomial of order k respectively. 



Proof. Consider the following cases: 
• Case 1: a is odd 
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Observe that 



JO 



2(_i){°-i)/2^ / sm{ex)e-' /^de 



2(-l) 



{a-l)/2 



d" 



dx"'' 



-Im 



2(-l)("-i)/2_^Im 



/>oo— la; 



by letting q = 9 — ix. We evaluate the integral 



J —1 



e-'^'/^dq 



by forming a closed anti-clockwise rectangular contour from —ix to R — ix to R to to —ix. 
The contour integral is zero by analyticity. Note that 



g2/2 

V 2 



'^dq 



The integral from R — ix to R tends to zero as -R — ?■ oo. By a change of variable t = q/i\^, 
the integral along the segment from to —ix is 

e-'fl'^dq = / e* 
JO 

= -iV2e-'V2^(^) 
where -D(-z) is the Dawson's integral. Therefore we have 

^e-V2,, = y| + ,V2e^V2^(^) 



and the result follows. 



Case 2: a. is non-integer 
Again we have 

/oo /•oo 
e-*^"-^ />rd0 = 2 / cos(0x)e-^ 
-oo Jo 

Using sum-of-angle formula and Q, we can write the classical parabolic cylinder function as 



./l COS — e^'/^ r e-*'/2t" cos(tx)dt + ^/^^sin — e^''/^ T e'^'/^t" sin(tx)dt 
V vr 2 7o V vr 2 7o 



TT 2 

Since the first term is even while the second term is odd in x, we have 

D^{x) + Da(-x) = 2a/ - cos —e"''!^ H e-*'/^^" cos{tx)dt 
V vr 2 Jo 



The result follows. 
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• Case 3: a is even 

Note that the integral now becomes analytic and the result is obvious by the discussion in 
Section 2.1. 

□ 

We are now ready to prove Theorem [2) 

Proof of Theorem\^ We follow the proof outline of Evans and Swartz (2000). First note that we 
can express the characteristic function 

where x s-nd ^ are as defined in (jl2p and ()13p . and R is the error term. Also let 



m = x{e) + -+m 



be the expansion of characteristic function but trimming the first term and the error term. 
Let 

I n, /9_1 I , 

{n\{e/^)f 



fc=0 



and qn{0) be the expansion of Wn{0) up to order 

f l/n"/2-iL(^) for 



non-even a 

i/2. 



l/n"/2-iCi(l^) for even 



a 



In other words, 



^ 1 . 1 I - r^^^iL. ill '^'"^i'/2-/'^'^ for non-even a 



4<j<afe=l r(a+l) n"/2-i 



for even a 



where are defined as in ([3]). Also let r„(0) = Wn{0) — qn{ff) be the error of this trimming. 

It is easy to verify through direct expansion that r„(0) is a sum consisting of terms written as 
Pk{9,n)n~^ , where k > a/2 — 1 and Pk{6,n) is a sum of terms in the form 

[ ' ) for non-even a 

A ( 1/2 \ 

|g| — -)'^ for even a 

for c G M, 6 > and d G N (and note that for A; = a/2 — 1 the terms in Pk{0, n) must have d = 0). 
We first show that 

-oo 
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is a good approximation to fg^/^{x) = e *^^+'^'^(^/v^)(i0/27r. After that we will verify using 
Lemma m that this approximation is the one shown in Theorem [2j Consider 



< 



oo 



de 

2^ 



2tt 



2tt 

+ 



27r 



(23) 



for any 5. Later we are going to choose the value of 5 to bound our error. 

We first consider the third integral in (j23p . By non-lattice assumption, given any 8, we can find 
< 1 such that \4>[y)\ < Tg for any y > 6. Hence 



de 
2iT - ^ 



n-l 



2^ 



o{n~P) 



for any p > 0. The last equation follows from our assumption that (p is integrable. 
Now consider the second integral in (|23|) . Note that for any /3 > 0, 

e~''/^9^d9 = o{n-n 
for any p > 0. Also, for any a > 2 and slowly varying function U, 



(24) 



let 



de < c„ 



e /2|^|"-%9/2^^ ^ q(^-p) (25) 



for any q > and Cq a positive constant depending on q. Again p > is arbitrary. By plugging in 
the expansion of together with observations (p^ and (p5]) . we conclude that the second integral 
in (j23|) is o{n~^) for any p > 0. 

Finally consider the first integral in (j23p 



/V(^)l- 



< 



2tx 



dB_ 
2^ 



+ 



( 



2/2 



^n(^)| 



Ci0 

2^ 



By our observation on r„(0) that it can be written as a sum of Pk{e,n)n~^ with conditions on 
and k, together with the fact that /_^^^e~^ ^"^{ei^de < qg for any 6 > and the slowly varying 
functions L and grow slower than any polynomial, we can verify that 



t 



2/2 



de 

2^ 



a/2-1 
,q/2-1 



for non-even a 
for even a 



Now consider the first part. We use the estimate 



j=0 



(3^ 



E 

j=0 



(3^ 



<i\a-P\ + 



17+1 



(7 + 1)! 



gmax{a,/3} 



(26) 
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where we use mean value theorem in the last inequality. Also, we can find 5 small enough so that 
for 9 G [—5,5], we have 



< - and \Xi9)\ < - 



which implies that for 6 £ [—5^/n,5y/n], we get 



Using the above estimates (j26|) and ([2 



< — and 
4 



.de 



nX 



- 4 



(27) 



< 



271 



nR 



e 



2/4 



(S\/n 



(La/2-lJ+l)! 



La/2-lJ+l 



Now since 



X{9) = o{9^+P) 

for ^ G [— 5-y/n, J-y/n] and some small p > 0, it is easy to see that the second integral in 



t 

-<5\/ri 



-6*2/4 



nX{e/^/n) 



(La/2-lJ+l)! 



La/2-lJ+l 



de 

2^ 



n 



a/2-1 



On the other hand, we have 



ar 1 



m 



o lei'^L 



for non-even a 
for even a 



(29) 



Suppose that a is non-even. For given e, we can choose 5 small enough such that for e G [—5, 5], 

R{e) < el^rL 
which implies that for e £ [—6^/n,6^/n\, 

\er 



nR 



and the first integral in (j28l 



ev4 



nR 



< e- 
n I n! 



de 



a/2-1 



L 



n 



1/2 



-02/4_ 



-L 



n 



1/2 \ d0 



n 



I 1^1 / 27r 



Note that by Lemma [2] we can pick M > and any Pi, P2 > such that 



^^1^^ ^ / C|6l|-^i for \e\ < 1 



Lin^^) - 1 C\e\P'^ for |6'| > 1 



for n > M. So by dominated convergence theorem and that 



.1/2. 



H 

L(nV2) 



as n — > oo for all fixed nonzero 9, we have 



n"/2-i" \^ 1^1 J 2tt 



f 

- 7_oo ' ' i^(nV2) 27r n"/2-i 

/oo 
-oo 



' ' 27r n"/2-i 



as n oo. Hence 



lim sup 



L(nV2)/n«/2-i 
for some constant C > 0. Since e is arbitrary, we get 



nR I ^ 



2^ = " 



'L(ni/2 



x) 



n' 



,a/2-l 



,a/2-l 



a is the same except L is replaced by that is also 

varying function. 



for any fixed x. The case for even 
:ying function. 

We can therefore conclude that 



/oo p 
-oo 



dd 



= < 



„a/2-l 



for non-even a 



y o I ^o,/2-i ) even a 

We now evaluate J^^ e''^"^-^^ qn{e) § . Note that 



f 

J — ( 



-iex-0^/2 



nJ72-i 

3<i<a 
j even 



/oo 
-oo 



idx-e'^/2 



V2 y 27r 



where Gj's are the ordinary Edgeworth coefficients. So it suffices to focus on the last term 



f 

J — ( 



-iex-eV2^^ 



9 \ dO 



IT 



nV2 ) 27r 

a+i)smif J-oo^ 1^1 ^il9r''2? n"/2-i lor non-even a 



2(-l)°/^ TOO i0^_e2/2m,a> /ni/2.de 1 for pvpti 

r(a-+-i) J-oo^ 1^1 '^L\-ifi\^)2i ;i^;72^ tor even a 
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By dominated convergence theorem and ([30]) . we have for odd a, 



vr 


r(a + 1) 


sin^f 


vr 




r(a + i) 


sin^f^ 


vr 





\9\ I 2vr n"/2-i 



-oo 
oo 



L(nV2)2vr n"/2-i 



r(a + 1) sin ^ 27r n"/2-i 

\/2 ^ X L(nV2x) 

by using Lemma H] in the last equahty. 

For non-integral and even a essentially the same argument follows, again noting that Cl(^) 
also a slowly varying function for proving the even a case. □ 



4 Beyond Central Limit Region 



Although this paper concerns expansions in the central limit region, we devote this section to an 
interesting connection of our extra term with the large deviations of regularly varying rv's. 
By local analysis (see for example. Bender and Orszag (1999)) 



dx" ^^2 ^/2xi+° 



as X ^ oo and for non-integral z^. 
So for odd and non-integral a, 



^"6 



27r g^V4|^|-^-l as Z 



as 2 — +00 
— oo 



-(l+a) 



as x — 7- OO. Note that the result for non- integral a follows from the observation that the negative 
side of Da in Ga{x) dominates and the reflection property of gamma function. Thus when a is 
non-even, we get 

L{n^/^x) L{n^/^x) 



Gaix) 



(31) 



1 j,l+a^a/2—l 
for X — > OO. 

On the other hand, Rozovskii (1989, 1993) proved a result for regularly varying rv on the whole 
axis: 



PiSn > X) 



+ nP{Xk > x) 



(l + o(l)) 



(32) 



uniformly for x > -y/n, where $(•) denotes the complementary normal distribution. In other words, 

P{Sn > X) 



lim sup 



(^{x/^/n) + nP{Xk > x) 



1 
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In particular, the result implies that for x > y^(a — 2)n log n, the large deviations asymptotic 
(the second term on the right hand side) dominates i.e. P{Sn > x) ~ nP{Xk > x). If we assume 
the existence of density, this can be written as 



fs 



x) ~ nf{x\/n) ■ y/n 



for x > Y^(a — 2) log n, which is exactly ()3ip . In other words, the non- analytic term matches exactly 
with the large deviations asymptotic. 

However, for even a our extra term resembles standard Edgeworth expansion, and there is no 
region on the real axis that this can coincide with the large deviation asymptotic. 

To rigorously show the validity of the connection would require knowledge about the second- 
order variation of the slowly varying function L i.e. how fast L{xt) / L{x) — )• 1 for given t. This will 
be a direction for future research, but in this paper we are satisfied that at least in special cases 
this connection is easily seen to be valid: 

Theorem 3. Assume {X^ '■ k > 1) are iid rv's with unit variance and symmetric Pareto density 
f{x) = af/{hf + Ixl""''"'^) where a is non-even, and af and hf are constants. Then we have the 
asymptotic 



1+ E 



3<j<a 
j even 



ni/2-i 



n 



a/2-1 



a/2-1 



n 



(33) 



uniformly over x < \/{a — 2) log n, where Gj (x) are defined as in Theorem [21 



In other words, the refined moderate deviations over the range x < (a — 2) log n is the same 
as the Edgeworth expansion together with our additional term. 

Proof. The proof is a simple adaptation from the proof of Theorem [2j Note that in that proof we 
bound uniformly the error of Edgeworth expansion over the real axis, but this error can be too 
coarse for x outside the central limit region. It suffices, therefore, to refine our error estimate so 
that it always has smaller order than our Edgeworth expansion over x < {a — 2) log n. In other 
words, we need to show that the error of our approximation using Edgeworth expansion is of order 
o{L{n^/'^x)/ {x^~^°'n°'/'^~^)) for any x < ^/{a^^2Ylogn. 

We will first prove the result for the odd a case. An easy examination of our proof of Theorem 
[2] would reveal that for Pareto rv, the only error that needs further refinement is that caused by 



-02/4 



nR 



dO 
2^ 



For this we need a finer estimate of R{9), which in turn needs a finer version of Case 1 in Lemma 
m Note that 



idx 



1) 



|l+o 



dx = 2a 



2ate 



2a/6l 



oo j.a— 1 



{cos{9x) 



I) 



-dx 



6/0^+" + 



du 



cos u~ 1 



du - 2a/ 6/^2+" 



(cos n — 1) 



7x2(6/^1+° + ni+" 



-du 
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The first part becomes —afirO, while the second part is 

(cos u — 1) 



n2(6/0i+" + 'ui+") 
(cos u — 1) 



du 
du + 



(cos u — 1) 



u- 



-du 



where e is a small number. Note that for < n < e we have (cosu — l)/ii^ bounded around — 1. 
Hence 



(cos n — 1) 



'u2(5^0l+" + ^l+a 



-du 



< c 



C 



1 



6/^1+" + ui 
1 



-du 



< 



9'^ Jo + 
C 



dx 



where C are positive constants. Obviously 

(cos u — 1) 



Therefore 



a/x"-i(e*''^ - 1) 



dx 



du = 0(1) 



-2a/7r6' + 0(6|2 



3 6/ + |x|i+" 

and hence -R(^) = O(|0|^~''"). Replacing (I29p with this estimate of -R(^) and applying the same 
argument as in the proof of Theorem 2 implies that 



<5v^ 

( 

-<5i/)i 



-^2/4 



d6l 
2^ 



O 



1 



a/2-1/2 



L(n^/2x) 



for X < y^(a — 2) log n. This leads to a uniform estimate of our error in the Edgeworth expansion 
that always has order less than Gaix)/n°'^'^~^ over x < ^/ {a — 2) log n. 

Now for non-integer a, when q = [a\ is even, Case 3 in Lemma [J can be written as 

a/a;«(e^^^ - 1) 



oo 



-dx = 2a/6'"-'? 



r-, ; — (cos u — l)du 



Now 



cos u — 1 



ni+"-'/(6/^i+" + ni+°) 



2a/ 



W/ 

°° COSU - 1 







u 



1+a-q 



du — bfi 



COS u — 1 



ui+"-'?(6/0i+° + 



cos n — 1 



-du + 



cos n — 1 



for small e > 0. Clearly the second integral is bounded. Consider the first integral 



cos u — 1 



ui+"-''(6^0i+" + 



du 



< C 



u 



Of 

C 



du 



1 



y 



l-(a~q) 







5 



O 
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where C is a positive constant, and a substitution u = 6y is used in the first equahty. Hence 
a^x'^^e^o- 1) ^^ ^ 2a.r-''r(-a + q) cos ( ^^n) + 0(9^) 

When q = [a J is odd, we have 



, ^ „^ dx = 2m.r-« 



&/^^+" + ni+' 



-du 



Now 







2iaf 



-du 



« 



smu 



du — bfi 



Sinn 



^i 



-du 



smu 



ui+°-9(6jei+" + ni' 



-du + 



smu 



for small e > 0. Clearly the second integral is bounded. Consider the first integral 

sinu 



u 



-du 



- Jo u^-i{hfe^+^ 



,0 u- ■i{bf9'+'^ + u^+^) 
C f-e 1 



-du 



-dy 



where C is a positive constant, and a substitution u = 6y \s used in the first equality. Hence 

°° arr'^e^^^ f -a + a \ 

^ = 2m;r-^r(-a + q) sin (^^^J + 



Therefore 

, , _ J O(|0| " when loj is even 

" \ 0(1^1""") when LOj is odd 

and the same argument as the odd a case holds. 



□ 



5 Non-Symmetric Density 

With the result for symmetric density in hand, we can generalize our results to the non-symmetric 
counterpart, namely Theorem [TJ The main idea is to write a non-symmetric density as 



where 



/(x) = r(x) -|- s(x) 



(34) 
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is a symmetric function, and 



s{x) = hf{x)-f{-x)] 



(35) 



is an odd function. The symmetric function r can be easily handled using the results in the last 
section, and the odd function s can also be dealt with by slight modification. Indeed we have the 
following proposition: 

Proposition 4. For rv following ^ with Lj^ and L_ satisfying Assumption\^ we can decompose 
its density as 

f{x) = r(x) + s(x) 

where r and s are defined by ^34\ ) and ^35\) respectively. Then r and s must be regularly varying 
and have the same order a = mm{(3,j), and we can write 



1 + kc 



1 + |x|i+" 



-CXD < X < OO 



(36) 



where is even and Lg is odd, and they are both slowly varying functions satisfying Assumption 
d Also define 



m. 



(r) 



/oo 
x^s{x)dx for < n < a — 1 
-oo 



as the "moments" of r and s. The cumulant generating function takes the form 

m = x{.o)+m+o{m) 

where 

l<j<a ■'' 

is the ordinary Taylor series expansion up to the largest cumulant, and 



m 



2tt 



r(o+l) sin(a7r) 



cos^L,(^)Tisin^fL,( 1 



for even a 

for odd a 

for non-integer a 



(37) 



is the non-analytic component. Here it and =F both refer to the cases when 6 > or 9 < 0. 



Proof. That r and s are regularly varying, have the same order and can be written as ()36p satisfying 
Assumption [1] is obvious. Now define (^r(^) (j)si^) the Fourier transforms, or "characteristic 
functions" of r and s. Since r is symmetric, by the same argument as Proposition [3l we get 

''^ ^ \h. ' ~ \ ^^irC.,(^)+o(lrC.„(^)) forevenc 

j even 

(38) 



28 



Now consider the odd function s. Note that the role played by odd and even a (or q) is reversed 
because of odd instead of even function. Also note that (f)g{0) = —(/)g{—6), so it suffices to prove the 
result for ^ > and the 6 < case follows. Following similar line of proof as Proposition [3l we get 

6(0)-!^ + / ^("+^rcos^ \^ns il) + o {\en, (i)) for non-odd a 

^'^ ^ ~ ^ ,t<P I rlriT l^rC. (^) + o (irC. (^)) for odd a ^''^ 

j odd 

for ^ > 0. Adding ()38p and (I39p . and noting that the representation can be carried over to cumulant 
generating function, we obtain the stated result. 

□ 

We also need a correspondence of Lemma H] for odd function, the proof of which is very similar 
and is thus omitted: 



Lemma 5. 



/oo 
-oo 



f 2^/2i-+i^Z)(^) for even a 

V2^e-'''/^i''Ha{x) for odd a (40) 

5 CSC ^e~^ /^[Da{x) — Da{—x)] for non-integer a 



where D{z), Hk{z) and Dy{z) are Dawson integral, Hermite polynomial of order k and classical 
parabolic cylinder function with parameter v respectively. Here 



1^1" /or > 
-l^l" fore<0 



The proof of Theorem [T] follows easily from Proposition H] and Lemma [5j 
Proof of Theorem [IJ By ([2]) we can write Lr and Ls as 

Lr{x) = UL+{\x\) + L^{\x\)] L,(x) = ±i[L+(|x|)-L_(|x|)] when /3 = 7 

L,(x) ~ iL+(|x|) L,(x)~±iL+(|x|) when/3<7 (41) 

Lr{x) ^L_{\x\) Ls{x) ^^L^{\x\) when /3 > 7 

where the it and =F signs correspond to x > and x < 0. We can proceed similarly as in the proof 
of Theorem [2] by substituting ([^T|) into Lr and Ls, and using Lemma O Note the use of ([7]) in 
obtaining the /3 < 7 and /3 > 7 cases. □ 



6 References 

1. Bazant, M. Z. (2006), 18.366 Random Walks and Diffusion, Lecture Notes, MIT OpenCourse- 



Ware, http://ocw.mit.edu/OcwWeb/Mathematics/18-366Fall-2006/LectureNotes/index.htm 



(See the notes and supplementary material for lecture 7.) 



29 



2. Bender, C. M., and Orszag, S. A. (1999), Advanced Mathematical Methods for Scientists and 
Engineers: Asymptotic Methods and Perturbation Theory, 2nd edition. Springer- Verlag. 

3. Bouchaud, J. P. and Georges, M. (1990), Anomalous diffusion in disordered media: Statistical 
mechanisms, models, and physical applications. Physics Reports 195 (4-5), 127-293. 

4. Bouchaud, J. P., and Potters, M. (2000), Theory of Financial Risks, Cambridge, UK: Cam- 
bridge University Press. 

5. Embrechts, P., Kluppelberg C, and Mikosch, T. (1997), Modelling Extremal Events for In- 
surance and Finance, Springer. 

6. Evans, M., and Swartz, T. (2000), Approximating Integrals via Monte Carlo and Deterministic 
Methods, Oxford Statistical Science Series. 

7. Field, C. A., and Ronchetti, E. (1990), Small Sample Asymptotics, Institute of Mathematical 
Statistics, Lecture Notes - Monograph Series. 

8. Hall, P. (1992), The Bootstrap and Edgeworth Expansion, Springer- Verlag. 

9. Hughes, B. (1996), Random Walks and Random Environments, Vol. 1, Oxford, UK: Clarendon 
Press. 

10. Metzler, R. and Klafter, J. (2000), The random walk's guide to anomalous diffusion: a frac- 
tional dynamics approach, Physics Reports 339, 1-77. 

11. Resnick, S. I. (1987), Extreme Values, Regular Variation, and Point Processes, Springer- 
Verlag. 

12. Rozovskii, L. V. (1989), Probabilities of large deviations of sums of independent random 
variables with common distribution function in the domain of attraction of the normal law. 
Theory Probab. Appl. 34(4), 625-644. 

13. Vinogradov, V. (1994), Refined Large Deviation Limit Theorems, Longman Scientific &; Tech- 
nical. 

14. Weiss, G. H. (1994), Aspects and Applications of the Random Walk , North-Holland. 



